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(57) A method of measuring internal information 
in a scattering medium, and an apparatus for 
the same of the present invention measure 
internal information in the scattering medium 
by measuring light diffused during propagation 
in the scattering medium while receiving the 
influence of the scattering and absorption, and 
performing arithmetic processing to the 
measured values. At this time, the three or more 
kinds of the detected signals (measured values) 
measured at three or more different kinds of 
distances between the light incident position 
and the photodetectfon point are processed by 
utilizing dependencies of the behavior of light 
diffused during propagation in the scattering 
medium and the resulting signal, that is, a 
photodetection signal on characteristics such 
as a scattering constituent, or an absorption 
constituent in the scattering medium and the 
concentration. The measured internal infor- 
mation is the absolute value of a scattering 
coefficient or a transport scattering coefficient 
of the scattering medium, and the information 
associated with a specified scattering con- 
stituent or a specified absorption constituent 
can be measured by processing the measured 
internal information. 
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BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to non-invasive measurement of internal information in a scattering medium 
by causing pulsed light, square wave light, or continuous light to be incident on a scattering medium such as 
a living body, and detecting light propagating in the scattering medium and more particularly, to a method for 
measuring internal information in a scattering medium and an apparatus for the same, which are capable of 
measuring an absolute value of an absorption coefficient and a transport scattering coefficient of the scattering 
medium, the concentration of a specific absorptive constituent in the scattering medium, its time change, its 
spatial distribution, and the others, and which are capable of improving measurement precision. 

Related Background Art 

is Demands for precise measurement of an absorption coefficient or a transport scattering coefficient in a 

scattering medium such as a living body or the precise measurement of the concentration of a specific ab- 
sorptive constituent are very strong, and there are several reports and attempts therefore. Main prior arts will 
be represented as references'- *> at the end of this section. Note that the number of each reference will be re- 
ferred to with a superscript*) hereinafter. 
20 in general, since light is scattered and/or absorbed at random in the scattering medium, light does not prop- 

agate straight For the scattering medium in which the absorption is zero, a total quantity of light never de- 
creases, but since light is scattered by scattering constituents at random, light propagates zigzag. In this case, 
a distance that light can propagate without the effect of the scattering is called a mean free path or a mean 
diffusion length, which is an inverse of a transport scattering coefficient u, s \ In the case of the living body spe- 
cimen, the mean free path is about 2 mm 1 ). In addition, the scattering medium comprises absorptive constit- 
uents other than the scattering constituents, so that the absorption occurs in accordance with a distance that 
light propagates while scattering, and the quantity of light is exponentially attenuated in accordance with the 
distance. 

As the prior arts in the field of the precise measurement of an absorption coefficient or a transport scat- 
tering coefficient of a scattering medium, one is to measure quantity of transmitted light or reflected light cor- 
responding to continuous light or pulsed light incidence, and one is to measure transmitted light or reflected 
light corresponding to pulsed light incidence and to analyze its waveform. The former utilizes a measurement 
of the absorbance in which Lambert-Beer Law is a basic principle, and there is one further utilizing a principle 
of dual-wavelength spectroscopy. Here, Lamber-Beer law is that the absorbance of the specimen is propor- 
tional to the product of the molar absorption coefficient, molar concentration, and a specimen thickness, and 
that a difference of the absorbances is proportional to the concentration difference where a thickness of the 
specimen is constant 

However, in the case of the absorbance optical density measurement in the scattering medium, a mean 
optical pathlength of light diffused during propagation between a light incident position and a photodetection 
point is varied depending on the absorption coefficient *i a of the scattering medium. Therefore, in the absor- 
bance measurement corresponding to the scattering medium in which the optical pathlength is constant, an 
absorption coefficient dependency of the mean optical pathlength is a big problem, which hinders the precise 
measurement of an absorption coefficient or the concentration of an absorptive constituent, or if the meas- 
urement is performed, due to large errors in measurement, it cannot practically be used. For example, accord- 
ing to "Japanese Patent Application No. Sho 62-248590", "Japanese Patent Application No. Sho 62-336197". 
and "Japanese Patent Application No. Hei 2-231 378". since a basic principle is to measure an optical coefficient 
in which an optical pathlength is assumed to be constant, errors in measurement caused by a change of the 
above-described optical pathlength cannot be avoided. 

There is another method utilizing a mean optical pathlength which is measured by another method in the 
case of an absorbance measurement, but hence a mean optical pathlength is varied depending on the absorp- 
tion coefficient, errors caused by approximating the average optical pathlength to be a constant value cannot 
be avoided. Further, there are a method of measuring an absorbance difference using pulsed light, a method 
further utilizing a principle of dual-wavelength spectroscopy, and a method of measuring an absorbance optical 
density using light having three or more kinds of wavelengths, but in either case, since a method of measuring 
an absorbance in which an optical pathlength in a scattering medium is assumed to be constant is applied, 
errors caused by the change of the optical pathlength occur in the all cases, so that the sufficient precise meas^ 
urement can not be performed. 

In the prior art not based on an absorbance measurement, a method, in which transmitted light or reflected 

2 - 



30 



35 



AO 



45 



50 



55 



EP0 627 620 A1 



light is measured by time-resolved measurement using pulsed light or modulated light to analyze waveform, 
has disadvantages that owing to the time-resolved measurement, a measurement method and an apparatus 
are very complicated and that the apparatus is expensive, as compared with a measurement method and an 
apparatus of the present invention in which quantity of light, that is, time integration of an optical signal is meas- 

5 ured. There are several attempts for measuring internal absorption information by measuring reflected light 
or transmitted light upon the incidence of pulsed light to the medium by the time-resolved measurement and 
analyzing its waveform 3 - 5 >. A measured optical signal has a long decay tail by the influence of scattering and 
absorption constituents. Patterson et al. assume a model of uniform scattering medium to analytically obtain 
the light signal output 3 ). A wave representing a time change in intensity of the optical output signal given by 

10 the formula defined by Patterson et al. matches a waveform obtained by an experiment using a uniform scat- 
tering medium. According to Patterson et al. and the results of the experiment by the inventors of the present 
application, the absorption coefficient of absorptive constituents in the scattering medium is given by a slope 
of waveform (differential coefficient) obtained when the optical signal is sufficiently attenuated, i.e., when a 
sufficiently long period has elapsed. However, because the optical signal at a location where an absorption 

is coefficient is obtained required to be sufficiently attenuated means that the signal is very feeble, a signal to 
noise ratio (S/N) of the signal to measured is decreased, and consequently errors in measurement are in- 
creased. Therefore, it is hard to use in practical application. There are several kinds of attempts other than 
the above, but neither of them gives the sufficient measurement precision. 

A method which utilizes the above-stated dual-wavelength spectroscopy measurement to a scattering me- 

20 dum has the following problems. In the absorbance measurement of a scattering medium is that an extinction 
coefficient is a sum of a transport scattering coefficient and an absorption coefficient on the definition. Since 
these parameters are treated at the same level, the absorption influence, e.g., an absorption coefficient cannot 
be measured while scattering and absorption influences cannot be separated. In general, this can be solved 
by using a principle of dual-wavelength spectroscopy measurement. In the dual- wavelength spectroscopy 

25 measurement an absorption coefficient is measured by using suitable two or more kinds of lights having a 
different absorption coefficient to an absorption constituent, and assuming that an scattering coefficient and 
a transport scattering coefficient of the two or more kinds of lights are the same, or if they are different, as- 
suming that the difference is very small, the scattering influence is eliminated from an absorbance optical den- 
sity difference corresponding to two or more kinds of lights to obtain the concentration of an absorption con- 

30 stituent or an absorption coefficient. This method has a disadvantage of the occurrence of errors caused by 
the assumption that scattering coefficients or transport scattering coefficients corresponding to different wa- 
velengths are equal. 
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SUMMARY OF THE INVENTION 

50 Measurement methods and measurement apparatus in the above-described field which have been devel- 

oped and reported have not satisfied the measurement precision yet. The inventors of the present application 
theoretically and experimentally analyzed and examined the behavior of light in a scattering medium which is 
a basis of the measurement in the above-described field, and found that the conventional analysis methods 
are insufficient to describe the complicated behavior of light in the scattering medium, and that the most of 

55 approximations used in the process of deriving a principle of the measurement of the conventional measure- 
ment methods and the measurement apparatus do not hold in the actual object to be measured and the meas- 
urement field, which results in large errors included in the measured value. It is an object of the present inven- 
tion to provide a new measurement method and a measurement apparatus which enhance the measurement 
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precision while the errors are drastically decreased based on the above analyses, examinations, and experi- 
ments. 

A method of measuring internal information in a scattering medium and an apparatus for the measuring 
the internal information such as an absolute value of an absorption coefficient or a transport scattering coef- 
ficient in the scattering medium utilizing new knowledge in which a change of optical pathlength is taken into 
consideration, that is, a plurality of simultaneous relations between photodetection signals (measured values) 
when light having a predetermined wavelength is diffused while propagation in the scattering medium and 
reaches at predetermined photodetection points and a plurality of known parameters. Further, information re- 
lated to absorption such as the concentration of a specified absorptive constituent is measured by performing 
the above measurement using lights having a plurality of wavelengths. 

A method of measuring internal information in a scattering medium of the present invention comprising 
(a) causing light having a predetermined wavelength to be incident on the scattering medium, (b) detecting 
light having the predetermined wavelength diffused during propagation in the scattering medium at detection 
points corresponding to three or more kinds of points where the distance between a light incidence position 
and a photodetection point is different to obtain the three or more kinds of the detected signals, and (c) ex- 
tracting internal information in the scattering medium by performing arithmetic processing to the detected sig- 
nals based on three or more kinds of simultaneous relations between a scattering characteristic and an ab- 
sorption characteristic on a diffusion-propagation path in said scattering medium when the light having the 
predetermined wavelength is diffused during propagation and reaches at said detection points, and the de- 
tected signals. 

Here, the light having the predetermined wavelength may be light having two or more kinds of wavelengths 
each having a different absorption coefficient to a specified absorptive constituent of said scattering medium' 
and the anthmetic processing may comprise a process of extracting primary internal information in the scat- 
tenng medium obtained by the measurement with one kind of the predetermined wavelength, and a process 
of extracting secondary internal information in the specified constituent of the scattering medium, utilizing the 
pnmary internal information corresponding to light having the two or more kinds of wavelengths. 

An apparatus for measuring internal information in a scattering medium of the present invention comprises 
(a) light-emitting means for emitting light having a predetermined wavelength, (b) light-incident means for caus- 
ing the light having the predetermined wavelength to be incident on the scattering medium, (c) photodetecting 
means for detecting light having the predetermined wavelength, diffused during propagation in the scattering 
medium at detection points corresponding to three or more kinds of points where the distance between a light 
incidence position and a photodetection point is different, and obtaining three or more kinds of the detected 
signals, and (d) arithmetic processing means for extracting internal information in the scattering medium by 
performmg arithmetic processing to the detected signals, based on three or more kinds of simultaneous rela- 
tions among a scattering characteristic and an absorption characteristic on a diffusion-propagation path in the 
scattenng medium when the light having the predetermined wavelength is diffused during propagation and 
reaches at the detection points, and the detected signals. 

Here, the light having the predetermined wavelength may be light having two or more kinds of wavelengths 
each having a different absorption coefficient to a specified absorption constituent of the scattering medium' 
and the arithmetic processing means may execute a process of extracting primary internal information in the 
scattenng medium obtained by the measurement with one kind of said predetermined wavelength, and a proc- 
ess of extracting secondary internal information in the specified constituent of the scattering medium, utilizing 
the pnmary internal information corresponding to light having the two or more kinds of wavelengths. 

A method of measuring internal information in a scattering medium, and an apparatus for the same of the 
present invention measure internal information in the scattering medium by measuring light diffused during 
propagation in the scattering medium while receiving the influence of the scattering and absorption, and per- 
forming arithmetic processing to the measured values. At this time, the three or more kinds of the detected 
signals (measured values) measured at three or more different kinds of distances between the light incident 
position and the photodetection point are processed by utilizing a characteristic that the behavior of light dif- 
fused dunng propagation in the scattering medium and the resulting signal, Le.. a photodetection signal depend 
on charactenstics and the concentration of a scattering constituent, or an absorption constituent in the scat- 
tenng medium. The measured internal information is the absolute value of a scattering coefficient or a transport 
scattering coefficient of the scattering medium, and the information associated with a specified scattering con- 
stituent or a specified absorption constituent can be measured by processing the measured internal informa- 
tion. If the measurement is performed at a plurality of places on an object to be measured, the spacial distrib- 
ution of the internal information can be measured, and if the measurement is performed at the same position 
and at different time, the time-change in the internal information can be measured. 

In particular, in the method and apparatus of the present invention, light having a predetermined wave- 
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length is caused to be incident on a scattering medium and light diffused during propagation in the scattering 
medium is detected at three or more kinds of distances between a light incident position and a photodetection 
point to obtain the three or more kinds of the detected signals (measured values). Next, the three or more kinds 
of detected signals (measured values) and the known parameters are processed to measure the internal in- 
5 formation in the scattering medium based on the three or more simultaneous relations of scattering and ab- 
sorption when the light having the predetermined wavelength reaches at the photodetection points, and the 
detected signals (measured values). The relation among the detected signals (measured values), the parame- 
ters, and the internal information in the scattering medium is disclosed in the present application for the first 
time. 

10 When the above measurement is performed using light having a plurality of different wavelengths, select- 

ing the predetermined wavelength to have the different absorption coefficients to the specified absorptive con- 
stituents in the scattering medium, the specified absorptive constituents in the scattering medium can be 
measured accurately. In the measurement of the present invention, since transport scattering coefficients and 
absorption coefficients at a plurality of wavelengths are measured at high precision, even in the case of spe- 

15 cimens or wavelengths in which the scattering coefficient is different at two wavelengths, the absorption coef- 
ficient of the scattering medium can be measured at high precision. The dual-wavelength spectroscopy meas- 
urement is utilized in the precise measurement e.g., the concentration of the specified material, the degree 
of oxygen saturation of hemoglobin and their spatial distribution and their time change. 

Therefore, a problem of errors in measurement caused by the change of optical pathlength and a problem 

20 of the scattering coefficient difference between two wavelengths in the dual-wavelength spectroscopy meas- 
urement are solved, and the measurement precision is drastically improved. Further, since the present inven- 
tion is not the time-resolved measurement method, but is a method of measuring the quantity of light, that is, 
the time integration value of a photodetection signal, the structure of the apparatus is very simple and as the 
utilization factor of light is improved, the measurement time becomes short, and the signal to noise ration (S/N 

25 ratio) is improved. Accordingly, the measurement or the monitoring of the concentration of oxygen in a brain 
of a human who is being operated on or exercising is made possible by the present invention. 

The present invention will become more fully understood from the detailed description given hereinbelow 
and the accompanying drawings which are given by way of illustration only, and thus are not to be considered 
as limiting the present invention. 

30 Further scope of applicability of the present invention will become apparent from the detailed description 

given hereinafter. However, it should be understood that the detailed description and specific examples, while 
indicating preferred embodiments of the invention, are given by way of illustration only, since various changes 
and modifications within the spirit and scope of the invention will become apparent to those skilled in the art 
form this detailed description. 

35 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a view showing behavior of light in a scattering medium. 

Fig. 2 is a view showing a way of obtaining quantity of detected light in an scattering medium. 
40 Fig. 3 is a graph showing relations between quantity of detected light l[p] and detection distance p. 

Fig. 4 is a view for explaining a principle of the present invention. 

Fig. 5 is a graph showing absorption spectra of respective biological materials. 

Fig. 6 is a view showing principle of measurement of a slab-like specimen. 

Fig. 7 is a view showing a configuration of an apparatus of the first embodiment. 
45 Figs. 8 and 9 are views respectively showing a configuration of switching and mixing lights having different 

wavelengths. 

Fig. 10 is a view showing a configuration of performing detection of lights having different wavelengths. 
Figs 11-14 are views respectively showing a method of causing light to be incident on a scattering medium. 
Figs. 15-17 are views respectively showing a method of receiving light. 
50 Figs. 18 and 19 are views respectively showing a method of amplifying a detected signal with low noise. 

Fig. 20 is a view showing a configuration of an apparatus of the second embodiment. 
Fig. 21 is a view showing a configuration of a main part of an apparatus of the third embodiment. 
Fig. 22 is a view showing a detailed configuration of a light incidence position and a photodetection point. 

55 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The embodiments of the present invention will be described in detail with reference to the accompanying 
drawings. The same parts are represented by the same reference numerals, and the repetitive description is 
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omitted. 



(1) Method of Measuring Interna] Information of Scattering Medium 

5 (1.1) Principle of Internal Information Measurement of Scattering medium 

Light is scattered, absorbed and diffused during propagation in a scattering medium such as a living body 
and some of light components emerges on a surface of the scattering medium. That is, light is scattered and 
absorbed in the scattering medium, but also can be transmitted therethrough. For example, Fig. 1 is a view 
w showing behavior of light in the scattering medium and shows a state in which light is incident on one point of 
the scattering medium to detect diffusing light (reflected light) at other point Ught incident on the scattering 
medium is scattered at random and spreads all over the scattering medium, but in Fig. 1 , only a track of photons 
detected by a photodetector is shown. In other words. Fig. 1 shows a track of photons which are utilized in 
practical measurement. 

is It is known that the behavior of light in the scattering medium is precisely described and analyzed in ac- 

cordance with a photon diffusion theory. According to the photon diffusion theory, a light pulse - ■ ient on 
the scattering medium is that its pulse width spreads as the light pulse is scattered, absorbed, anc diffused 
during propagation in the scattering medium. Fig. 2 shows installation of a virtual light source to obtain quantity 
of light in the scattering medium. As shown in Fig. 2. when light is incident on a point (p=0, z=0) of a scattering 
20 medium surface, light to be detected at a point P (p, z), i.e., an optical signal is derived from a photon diffusion 
equation. In this case, since the photon diffusion is not present on the surface of the scattering medium and 
outside thereof, a boundary condition is required to be set to satisfy this requirement Patterson et at develop 
a method which satisfies the boundary condition by assuming a point light source of negative polarity and re- 
port that the theoretical values obtained from this method and the experimental values are well matched. In 
25 Fig. 2, the point light source of negative polarity is set at a o mark (p=0, z=-Zo). 

On the other hand, the behavior of each photon diffused during propagation in the scattering medium can 
be analyzed, experimented, and examined in accordance with Monte Carlo calculation with use of a computer. 
In addition, the experiment can be conducted with a physical model of a scattering medium or a living body 
specimen. From the recent knowledge, the results derived from the photon diffusion theory, the results of 
Monte Carlo calculation, and the results of experiment with a sample are well matched. Therefore, it can be 
considered that the behavior of light in the scattering medium is sufficiently described by the photon diffusion 
equation. The various analyses, experiments, and examinations conducted by the present inventors and the 
recent accurate analyses and experiments of the behavior of light in the scattering medium confirm that the 
behavior of light in the scattering medium can precisely be described by the photon diffusion equation. 

The photon diffusion equation which describes-the behavior of light in the scattering medium as described 
above is. for example, expressed by a following equation using a photon f luence rate ij>. 

(1/c) <dcj>(r,t)/dt) - DV2<><r,t) + M>(r,t) = S(r,t) (1,1) 

where 

the Photon fluence rate at position r, at time t [photon/mm2«sec] (r is vector) 
40 D: the diffusion coefficient [mm] 

u** the absorption coefficient [mm- 1 ] 

c: tne speed [mm/sec] of light in a scattering medium (given by a refractive index) 

S(r,t): light source [photon/mm^sec]. 

Since an inpulsed light source can be expressed by a delta function, a light inpulse incident on an origin 
45 (r=0) at t=0 can be expressed by the following equation. 

S(r, t) = 8(r, t) = 8(0, 0) = 5(0)#5(0) (1 .2) 
Consequently, the photon diffusion equation corresponding to inpulsed light incidence is as follows: 
(1/c) (&Kr.t)/at) - DV2<t»(r.t) + M>(r,t) = 6(0.0) (1.3) 
where 6(0,0) is a light inpulse incident on an origin (r=0) at t=0. 
50 The various optical constants on a scattering medium are, for 

u. s - the scattering coefficient [mirr 1 ] 

u s ': the transport scattering coefficient [mnr 1 ] 
Ht r : the transport attenuation coefficient [mnrr 1 ] 

\ieff'. the effective attenuation coefficient [mnr 1 ] 
55 9- the mean cosine 0 of the scattering angle 0, 

D = [3(u a * u s ')J- 1 = (3u, r )- 1 (1.4a) 
Ms' = (1-g)u s (1.4b) 

Utr = m + US = Ma + (1 - g)u S (1.4c) 
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Heff = [3u a (u a + VsW 2 = {3u a [u a ♦ (1 - ^Hs]}" 2 . (1 -4d) 

In the case that a spot-like light pulse is incident on a semi-infinite scattering medium, a boundary condition 
is satisfied by providing a point light source of negative polarity at a location (p=0, z=-Zo) in accordance with 
Patterson et al. In general, z© is about 1/u. s \ but strictly speaking, it varies depending on a light incidence meth- 
5 od or characteristics of scattering constituents in a scattering medium. This is confirmed by the inventors of 
the present application with use of Monte Carlo calculation. 

The photon diffusion equation (1.3) is solved under the above boundary condition, and an optical signal 
R(p, 0, t) [photon/mm^sec] at an arbitrary location (p, 0) of the surface of the scattering medium is derived. 
R(p, 0, t) = (4*D) - *2Zot - «2 x exp[ - (zo 2 + p 2 )/(4cDt)] x exp( - cu^t) (1 .5) 
10 Since the detected quantity of light l[p] is a result of the time integration of R(p, 0, t), it is expressed by 

P 

l[p] = JR( o , 0, t)dt (1.6) 
and this equation can be solved with use of integral formulas to obtain 
f 5 I [p] = zolu^zo 2 + p2)i/2 + 1]/[2rc(Zo 2 + p 2 )* 2 } x exp[ - H^zo 2 * P 2 ) 1/2 1 (1 .7) 

Taking the natural logarithm of l[p], 

lnl[p] = ln[Zo/(27i)] - (3/2)ln(Zo 2 + p 2 ) * In^z** + p 2 ) 1/2 + 1 ] - HefKzo 2 + p 2 ) 1/2 (1 .8) 
Therefore, when Zo is a known value, a measurement value l[p] and a known value p are substituted into equa- 
tion (1.7) or (1.8) to derive u^f in a scattering medium. Fig. 3 shows a relation between the measurement value 
20 l[p] and a detection distance p. Here, ^ s ,= 0-45 and Zo=1/0.45, and u. a is shown as a parameter. 

However, z<> in the above-stated equations (1 ,7) and (1 ,8) is actually Zq * *i s ' • Tnis 030 easily be understood 
with a next example. In the general measurement system, light is incident on a scattering medium from air, 
and light diffused during propagation in the scattering medium is detected in air. In this case, the above- 
described boundary condition, i.e., Zo is a value which depends on the refractive index difference between 
25 the air and the scattering medium. According to the accurate analyses and Monte Carlo calculation by the in- 
ventors of the present invention, for Zo=a/(ii a +tis'), it is confirmed that in the general measurement state, the 
value of a is within the range of 0.3 < a < 3. 

Further, when light is incident on the scattering medium or when light emerges from the scattering medium, 
if there are some light attenuation, for example, absorption by a surface colored layer, the optical signal R ex- 
30 pressed by equation (1.5) and the detected quantity of light expressed by equation (1 .7) or equation (1 .8) are 
decreased in proportion to the attenuation. In general, the emerged light is detected by a photodetector. There- 
fore, the detected signal (measured value) is a value related to the quantum efficiency of the photodetector. 
However, the quantum efficiency can be considered together with the above light attenuation factor. 

A measurement method which can obtain the high measurement precision even though there is the light 
35 attenuation will be explained below. Assuming the attenuation factor of the light attenuation is r|, the meas- 
urement value Q is derived from equation (1.8) to obtain, 

InQfp] = ln(W(2n)] - (3/2)ln(z 0 2 + p 2 ) + H\^f (Zo 2 + P 2 ) 1 ' 2 + 1] - m^o 2 + P 2 ) 1 * (1.9) 

for 

^=3/013 + ^)- (1.10) 

40 Here, "a" is a constant which is determined experimentally or based on the experience, typically 0.3 < a < 3. 

The unknown values in equation (1.9) are apparently u. a , u. s ' . and r\ t consideringa is a constant Accordingly, 
the photodetection distance p which is a controllable known value is set to three or more different values, and 
three or more measurement values Q are obtained, and based on three or more simultaneous equations as a 
function of three or more kinds of Q, the three unknown values, Le., the light attenuation factor t|, the absorption 

45 coefficient \i a% and the transport scattering coefficient u. s * can be calculated. In the present invention, the 
above-described basic knowledge is developed in the following application. 

Fig. 4 is a view for explaining a method of measuring internal information in a scattering medium of the 
present invention. Lights are detected at three kinds of detection distances to obtain the measurement value 
Q[p]. Note that a distance between the light incidence position and the photodetection point is simply called a 

so detection distance hereinafter. Three kinds of the measurement values at the detection distances p 1f p 2 , and 
p 3 are assumed to be expressed as 

Q[pil • tne measurement value at the detection distance p n 

QfpJ : the measurement value at the detection distance p 2 

Q[p 3 ] : the measurement value at the detection distance p 3 . 

55 Then, the three measurement values at the three kinds of the detection distances satisfy the following three 
simultaneous equations. 

lnQ[pi] = lnfozo/(2jr)] - (3/2)ln(zo 2 + P i 2 ) + In^Zo 2 + Pi 2 ) 172 + 1] - ^K^o 2 + Pi 2 ) 1/2 0-Ha) 
InQtpJ = lnfozo/(2701 - (3/2)ln(zo 2 + p 2 2 ) + ^[u^o 2 + P 2 2 ) 172 + 1] - Meff(zo 2 + P2 2 ) 1/2 (1.11b) 
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InQIpsJ = ln[W(2n)] - (3/2)ln(z<,2 + p32) + InbvK^ + pfyn + ^ . + p;)2)1/2 (1 11c) 

The three equations constituting the simultaneous equations (1.11a) to (1.11c) (hereinafter also generally 
called (1. 11)) are independent from each other, and the unknown values are p,, n s \ and n . Therefore the three 
unknown values m, p s ' and n . and p*, or Zo if desired can be obtained by using the three measurement values 
and p„ p2, and p 3 which are the known values or measured by another method. Note that in general „ and z„ 
are not needed to be obtained. To obtain p, and p s ', or p^, the simultaneous equations may be of any forms 
as long as the equations are independent from each other and derived from equation (1.6) or (1 7) When n 
can be treated as a known value, p, and ps\ and p^ can be obtained by solving the simultaneous equations 
of twoequations arbitrarily from equations (1.11)orthe simultaneous equations of two independent equations 
derived from equation (1 .6) or (1.7). The computation to solve such simultaneous equations can be performed 
at high speed using a computer. It is apparent that the computation precision is improved as the number of 
independent measurement values Q is large. 

In the above-described measurement, using lights having different wavelengths X, and X 2 . the absorption 
coefficients ^ and p^ with respect to each wavelength and the transport scattering coefficients il,,' and 
P32 with respect to each wavelength can be obtained. Therefore, for example, the degree of oxygen saturation 
of hemoglobin is calculated from a value of p^/p^. In such dual-wavelength spectroscopy measurement shown 
in Fig. 5. wavelength dependency of the absorption coefficient of the absorptive constituents contained in the 
scattering medium is utilized. For example, in the measurement of oxyhemoglobin and reduced hemoglobin 
oroxymyoglobin and reduced myoglobin, the spectra of which are shown in Fig. 5, light having a wavelength 
in which the absorption coefficient difference between oxidation and deoxidation is large, that is, light havinq 
700 nm to 1 .2 pm wavelength is frequently used. 

The significant point is that in the measurement method and apparatus according to the present invention 
he absolute values of the transport scattering coefficients m,' and p^ can be obtained. Therefore, even* 
though the assumption on the conventional dual-wavelength spectroscopy, that is, the assumption that regard- 
ess of the sameness or difference of the scattering coefficients of the scattering constituents with respect to 
lights having two different wavelengths, the scattering coefficients are extremely small, does not hold the 
measurement method and apparatus of the present invention can perform the measurement with high preci- 
sion. As descnbed above, the measurement method and apparatus according to the present invention can 
largely .mprove the measurement precision by expanding the application of the conventional dual-wavelength 
spectroscopy measurement. 

The present invention utilizes the integral value of the detected optical signal, so that the incident light pulse 
waveform of which is of any forms can be utilized without any problems as long as the integral time can be 
specified. For example, the square wave or continuous light is apparently suitable. If the measurement is per- 
formed with lights having three or more kinds of wavelengths or at four or more kinds of detection distances 
improvement in the measurement precision and/or measurement of specimens which have background ab^ 
sorption are possible. In Fig. 4, lights are detected at the three different points for one light incidence position 
but a method in which light is incident on three different positions and light is detected at one photodetection 
point may be used. In a word, three or more kinds of the detection distances should be used. 

The case that the object to be measured is a semi-infinite scattering medium has been explained however 
in practice, a finite scattering medium is mostly measured. In this case, the boundary condition should be sat- 
isfied at the surface and the outside of the scattering medium shown in Fig. 4, and the condition of photo dif- 
fusion should hold in most of lights which are diffused during propagation in the scattering medium. This con- 
dition holds when the scattering medium is considered to be sufficiently lame as compared with Pl p, and 
P3 . For example, the region in the wavy lines in Fig. 4 shows a region through which most of lights pass' The 
presence of such a region is apparent from a spindle-shaped beam in the scattering medium reported by Sevick 
etal. or the results of Monte Carlo calculation. 

When the thickness (z direction) of the scattering medium is not considered to be sufficiently large as com- 
pared with the photodetection distances Pl , p 2 and p 3 , second point light sources of negative polarity and pos- 
itive polarity which satisfy the boundary condition may be installed on the back surface (surface in opposition 
to light incident plane) of the scattering medium. In this case, theoretically, in order to compensate the influence 
of the second point light sources of negative polarity and positive polarity, third point light sources of negative 
polarity and positive polarity are required and further in order to compensate the influence of the third point 
hght sources, fourth point light sources of negative polarity and positive polarity are required. However since 
the influence of the installed point light sources to the emitted light rapidly decreases in order, the influence 
can be approximated by finite point light sources. In this case, the above-stated equations (1 5) to (1 9) are 
required to be modified by adding terms related to the influence of the additional point light sources of negative 
polarity and positive polarity. However, since this modification is only to increase the number of terms not to 
change the number of the unknown values and the independence of every equation, the same relation'as de- 
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scribed above can hold. 

The measurement method and apparatus of the present invention can also be applied to a slab-like spe- 
cimen, and its state is shown in Fig. 6. Light is incident on three points in Fig. 6, but light may be incident on 
one point and lights are detected at three points of different detection distances. In the case of the slab-like 

5 specimen, assuming the detection distance is d, a result derived from the photon diffusion equation (1 .3) which 
is equivalent to the above-stated equation (1 .7) is, for example, 

l[d] = exp( -dm ff )/(2^d3) x {(1 +3zo/d)[u^<d - Zo)+ 1] (d -Zo)exp(zomf f ) - (1- 3zo/d) 
hWd + Zo) + 1](d ♦ Zo)exp( - zou^)} (1.12) 
The equation which is equivalent to equation (1.9) becomes 

10 Q[d] = ril[d]. (1.13) 

However, these equations are modification forms of equations (1 .7) and (1.9), that is, the number of unknown 
values and the number of known values are the same, and the parameters are also the same, so that similar 
to the above case, the arithmetic operation is conducted and the same results as the above case can be ob- 
tained. Therefore, in the same way as described above, the three unknown values \i at u. s * and r\, and also m ff 

15 if needed can be obtained with use of the three measured values and the known values p 1f p 2 , and p 3 . Note 
that in general, r| and Zo are not needed to be obtained. When tj is the known value, u. a , ji s , and can be 
obtained from the simultaneous equations of arbitrary two independent equations. Further, the internal infor- 
mation relating to absorptive constituents such as the degree of oxygen saturation of hemoglobin can be ob- 
tained. In the same way as described above, the computation to solve the simultaneous equations is performed 

20 at high speed using a computer. 

Further, the internal information in the scattering medium as described above is a mean value of informa- 
tion along with the optical path of light which is diffused during propagation to the photodetection point from 
the light incidence position. Therefore, when the internal information is coarsely diffused enough as compared 
with the distance between the light incidence position and the photodetection point as shown in Fig. 4, or as 

25 compared with the interval of the light incidence positions, the spatial distribution of various kinds of the meas- 
urement values can be imaged, that is, a simple imaging technique. It is apparent that if the above-described 
measurement is performed at different time, the time change in the internal information can be measured, and 
it is applied to the measurement or monitoring for the amount of oxygen in brain or others. In these cases, the 
imaging and the computation can be performed at high speed using a computer unit which comprises a mem- 

30 ory, a display or others. 

Errors in measurement in both cases of a method of the present invention and a conventional absorption 
method in which the optical pathlength is assumed to be a constant will be explained to show one example of 
the effectiveness of the present invention. Considering the scattering medium which occupies infinite spaces 
as shown in Fig. 4, the optical signal to be detected is expressed by equation (1.5). The average optical path- 

35 length <L>, that is, the centroide of time-resolved waveform is considered for such a signal to obtain, 
<L> = [cj tRj(p, 0, t)dt]/tfR(p, 0, st)dt] = (3/2) (n a + » s ') (zo 2 + p 2 )/[1 + (z<> 2 + P 2 ) 1 *^,]. (1.14) 
Assuming that fi s *=0-45[mm- 1 ], 2o=1W as typical values for the living body specimen, and that p=50[mm] 
as the measurement condition, and that the average optical pathlengths are <L(^ a =0.01)> and <L(n a =0.02)> 
when n a is changed from 0.01 to 0.02, respectively, and that the optical pathlength is constant, error A in the 

40 average optical pathlength is given by 

A = {<L(n a = 0.01)> - <L(u. a = 0.02)>}/<L(u. a = 0.02)> x 100 = 14 (%). (1.15) 
Consequently, the absorbance optical density obtained under the assumption of the optical pathlength be the 
constant has about 14% error. To the contrary, the method and apparatus of the present invention perform 
the arithmetic processing utilizing the relation in which the change of the optical pathlengths is taken into con- 

45 sideration, so that the errors as described above can be eliminated. The above description is the same in the 
case of the measurement of the slab-like scattering medium, and there is small difference in the percentage 
of errors, but the effectiveness of the present invention is apparent 

(1.2) Measurement Method of Absorption Information 

50 

A method in which information relating to absorptive constituents is processed and extracted with use of 
the absorption coefficient obtained as described above will be explained below. 

(i) Measurement Method of Degree of Oxygen Saturation of Hemoglobin 

55 

Main absorptive constituents in a mammalian brain are water, cytochrome, oxyhemoglobin, and reduced 
hemoglobin. Absorption of water and cytochrome in a near-infrared range is as small as negligible with respect 
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to an oxyhemoglobin and a reduced hemoglobin. The oxyhemoglobin and the reduced hemoglobin have dif- 
ferent absorption spectra, as shown in Fig. 5. The skull is regarded as a scattering medium with respect to 
near-infrared rays. ~ 

. . 71 ? two ki " ds of "9 hts havin 9 wavelengths X, and X 2 are considered. Absorption coefficients correspond- 
ing to the wavelengths X, and X 2 can be obtained in accortance with the Lambert-Beer law as follows: 

fei = SHb,i[Hb] + eHto.iDHbO] (2.1a) 
where " ^ = CHbJHb] + eHbo^IHbO] (2.1b) 

8hw: the molar absorption coefficient [mm-i. K/M] of the reduced hemoglobin at the wavelength X, 

eHbo.,: the molar absorption coefficient [mm-i- IYH] of the oxyhemoglobin at the wavelength X, 

W the molar absorption coefficient [mm-i- M-i] of the reduced hemoglobin at the wavelength X 2 

he molar absorption coefficient [mm-i- M-i] of the oxyhemoglobin at the wavelength X, 
[HbJ : the molar concentration [M] of the reduced hemoglobin 

[HbO] : the molar concentration [M] of the oxyhemoglobin 
Since the degree Y of oxygen saturation is given by 

Y = [HbO]/([Hb] + [HbO]), (2.2) 

the following equation is obtained. 

Th~^ . . Ka1 ^ 32 = [Et * 1 + Y <«W>,1 - Wl/ISHM + Y&HbO* - (2.3) 

h .n2 Jnfn ^ r ( ; oncentration <* the ^"ced hemoglobin [Hb], the molar concentration of the oxy- 
hemoglobin (HbO] the molar concentration of the hemoglobin [Hb] + [HbO] and the degree Y of oxygen satur- 
ation are calculated using m, and Ma2 , or m,/^ obtained from the above measurement, and the known para- 
meters £Hbo,t, e hm and e^. If a wavelength (= 800 nm, isosbestic wavelength) which has the same 
absorption level for both an oxyhemoglobin and a reduced hemoglobin is used, the above equation can be 
made simpler. 

(ii) Presence of Background Absorption 

| e v>u a J^h' OUnd a , bS0 ^ ti0 " may n0t te ne 9 |ected in a livi "9 body. In this case, if the background absorption 

L^h. ih « T T Uk2 ' Xi ^ d6f ined 38 32 ' and a * "*pec«ve.y. the following equations are 

established in accordance with the Lambert-Beer law. 

Mai = e Hb.i[Hb] + e Hbo1 [HbO] + a, (2.4a) 
Ha2 = SHbJHb] + e H60>2 [HbO] + a 2 (2.4b) 
H a3 = em>.3[Hb] + e H& o.3[HbO] + a 3 (2.4c) 
The above equations are rearranged to obtain the following equation 
<M.i - M*)/(^ - Haa) = [<*«,, -.<*«) * Y( £Hbo>1 - EHb0i2 . £HD1 + ^ + K(ai . a2)] + [(£Mb3 . 

for 8 + 0,3 " £Ht, °' 2 " 8mm + 6Hb -^ * K(a * " a 2>l ( 2 - 5 > 

K = 1/([Hb] + [HbO]). (2.6) 

wStZn !£!!X££**1 T r itab,y se,ected so that ****** in 1,16 same way as described a ^. 

Hb], [HbO]. [Hb]) + [HbO], and the degree Y of saturation can be obtained by obtaining p a1 , ^ and Ma3 using 
the measurement values at the lights having three kinds of wavelengths and the known parameters, and s2 
strtut.ng j^. j^, Ma3 and the known parameters into equation (2.4) or equation (2.5). Note that the above con- 
livin^body 2 ^ 33 ^ a ° hieVed by Pr ° Per ' y se,ectin 9 a wavelength for the object to be measured such as a 

(2) Apparatus for Measuring Internal Information in Scattering Medium 
(2.1) First Embodiment of the Apparatus 

The first embodiment of the apparatus carrying out the method of measuring internal information in a scat- 
tering medium according to the present invention is shown in Fig. 7. A light source 2 using a laser diode or the 
like generates ights having certain wavelengths X, and X 2 . The wavelength of light from the light source must 
be appropnately selected ,n accordance with an object to be measured. Light having a wavelength of 700 [nm] 
o LTin IZ 9 Z * V P ^ ra K asSOciation with absorption of a hemoglobin or the like. Since the oxyhemo- 
globin and the reduced hemoglobin have different absorption coefficients, as shown in Fig. 5, the oxyhemo- 
globin ,s distinguished from the reduced hemoglobin by selecting an appropriate wavelength and thenmeas- 

^ 1 'f SOUrCeS SUCh 38 3 Mght emittin9 diode - a laser diode - an HeNe laser or others can 
be used. Further, the light source which generates pulsed light or square wave light as described above may 
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be used. 

A wavelength of light from the light source 2 is selected by a wavelength selector 4 and the light is incident 
on a surface of a scattering medium 20 which is an object to be measured through a light guide 6. In this case, 
there is another method utilizing a condenser lens or a pinhole, which will be described later. Since the average 

5 diffusion length is approximately 2 [mm] in the scattering medium such as the living body specimen as descri- 
bed above, the incident light is scattered before it propagates about 2 [mm], and loses the original direction. 
Therefore, the influence of the average diffusion length in a few cm scattering medium can be neglected, so 
that spot-like light must be caused to be incident on the scattering medium. Thick beam-like light may be 
caused to be incident on the scattering medium. This case is considered as a plurality of spot-like light sources 

10 in a line. 

A space between the light guide 6 and the object 20 to be measured is very small in the embodiment of 
Fig. 7. However, this can be made large, and a liquid medium or a jelly-like object (hereinafter called an inter- 
face material) having substantially the same refractive index and the scattering coefficient as the scattering 
medium 20 which is an object to be measured may be filled in this space. Since the light is diffused during 

15 propagation in the interface material and then incident on the object to be measured, no problems arise. If the 
surface reflection on the scattering medium is a problem, the influence of the surface reflection or others can 
be reduced by properly selecting the interface material. 

The light diffused during propagation in the scattering medium 20 is received by light guides 8, 9 and 10 
provided on the position spaced apart from the light incident position by distances p 1t pa, and p 3 , respectively. 

20 The interface material may also be used for the same reasons as described above. 

The first photodetector 12, the second photodetector 13, and the third photodetector 14 convert the re- 
ceived optical signals into the electric signals, amplify the signals and output the detected signals a, b, and c, 
respectively. Various kinds of photod electors such as a photomultiplier, a phototube, a photodiode, an ava- 
lanche photodiode, or a PIN photodiode can be used as the photodetectors 12, 13, and 14. Any photodetector 

25 can be selected as long as it has spectral sensitivity to detect light having the specified wavelength. If an optical 
signal is very feeble, a high-sensitivity and high-gain photodetector is used. The other parts of the photode- 
tector other than the light incident plane may preferably have the structure which absorbs or shields light. If 
light diffused during propagation in the scattering medium has light components having a plurality of wave- 
lengths, a wavelength selection filter (not shown) is inserted between the photodetectors 12, 13, and 14 and 

30 the scattering medium 20. 

An arithmetic processing unit 16 processes the detected signals a, b and c from the first detector 12, the 
second detector 13, and the third detector 14, and converts the signals a, b, and c into the measurement values 
d, e and f in proportion to the detected quantity of light, respectively. In particular, the arithmetic processing 
unit 16 integrates the detected signals a, b, and c in time by utilizing a signal g which synchronizes with the 

35 light generation of the light source 1 2 to calculate the measurement values d, e, and f in proportion to the de- 
tected quantity of light. Note that when the pulsed light is utilized, the synchronized signal g may be omitted. 
Next, the measurement values d, e and f, and the known parameters p 1f p 2 , and p 3 which are set or measured 
by another method are processed based on the three simultaneous equations corresponding to the wave- 
lengths X. 1f and X 2 to calculate the absorption coefficients and the transport scattering coefficients correspond- 

40 ing to the wavelengths X n and X 2 , i.e. , ji a1 and u. s1 ', and u.^ and m 2 '- This sort of arithmetic processing can 
be operated at high speed by a micro computer installed in arithmetic processing means. If desired, the con- 
centration of the specified absorptive constituents is calculated by using the obtained internal information, i.e., 
u. a1 and iia2. If background absorption is present, the lights having three or more kinds of wavelengths as de- 
scribed above are used. 

45 In the above case, there are a method using time-divided light having a different wavelength emitted, and 

a method using light having different wavelengths at the same time. As means for selecting a wavelength, there 
is a method (Fig. 8) utilizing a light beam selector using a mirror, a wavelength selector using a filter, or a light 
selector using an optical switch. In addition, a method of selecting a wavelength by a wavelength selection 
filter provided in front of a light incident position with lights having different wavelengths being formed into a 

50 coaxial beam, a method of causing lights having different wavelengths to be incident on the scattering medium 
in parallel (Fig. 9) and selecting a wavelength by a wavelength selection filter provided in front of the photo- 
detector (Fig. 10), and further a method of detecting lights having two kinds of wavelengths at a respective 
detection point using six photodetectors, and a method of properly combining these methods, are available. 
A means for causing the light to be incident on a scattering medium such as a living body may be constituted 

55 by a method utilizing a condenser lens (Fig. 11), an optical fiber (Fig. 12), or a pinhole (Fig. 13), or a light incident 
method using a gastro camera (Fig. 1 4) or the like, other than the method using light guides as shown in Fig. 
7. A means for receiving and detecting the light propagating through the scattering medium may be constituted 
by a direct photodetection method (Fig. 15), a method of detecting the light through an optical fiber (Fig. 16), 
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a method of detecting light through a lens (Figs. 17), or the like. 

If the signal detected by the photodetector is required to be amplified with low noise, a narrow-band am- 
plifier (Fig. 18), a lock-in amplifier (Fig. 19) or the like can be used. When the lock-in amplifier is used as a 
reference signal, the aforementioned synchronizing signal g is utilized. This method is effective for the meas- 
urement of a high dynamic range using square wave light or pulsed light If the absorptive layer is present on 
the surface of the scattering medium 20 such as a living body specimen, the measurement values d, e and f 
can be compensated using the absorption values measured at the position of the detection distances o, o, 
and p3. 

The light incident position and the photodetection points on the scattering medium are concurrently scan- 
ned (not shown), and internal information of every location of the scattering medium is obtained and accumu- 
lated in a frame memory (not shown), and then read by a TV system, whereby the image showing distribution 
of interna information can be obtained. If the measurement is performed at different time, the time-change 
in internal information can be measured. The arithmetic processing unit 16 has a function of storing thus ob- 
tained internal information, and a display recording means 18 shown in Fig. 7 displays and/or records the 
course of processing and these results. The arithmetic processing can be operated at high speed by a computer 
unit comprising a memory, a display or others. 

Moreover, the same measurement can be made to the slab-like specimen as described above In this case 
the light incident position and photodetection points are required to be provided on a surface opposing to the 
scattenng medium and set on the different detection distances so as to detect the transmitted light diffused 
during propagation as shown in Fig. 6. 

(22) Second Embodiment of the Apparatus 

The second embodiment of the apparatus carrying out the method of measuring internal information in a 
scattenng medium according to the present invention is shown in Fig. 20. This apparatus is for the measure- 
ment of the optical constants in the scattering medium, and the specimen 20 put in a container 24 is a scattering 
medium and its absorption coefficient and its transport scattering coefficient are unknown. The apparatus of 
this embodiment uses light having a predetermined wavelength, and its operation is substantially the same 
as the first embodiment except for use of dual-wavelengths. In Fig. 20, the parts which has the same function 
as descnbed in Fig. 7 are represented by the same reference numerals. 

Light having a predetermined wavelength from a light source 2 is incident on a liquid specimen 20 which 
•s put in the container 24 through a light guide 6. Light propagated during propagation in the specimen which 
is an object to be measured is received by light guides 8, 9 and 10 placed at positions spaced apart from the 
light incident position by distances p„ p 2 and p 3 , converted into the detected signals a, b and c by a photode- 
tector (not shown) provided in a detecting unit 22, and provided to an arithmetic processing unit 16 in the next 
stage. The arithmetic processing unit 16 processes the detected signals a. b and c, and converts the signals 
into the measurement values d. e and f in proportion to the detected quantity of light, respectively. Next the 
measurement values d, e and f, and known parameters Pl , p 2 and ^ measured or set by another method are 
processed based on the simultaneous equations as described above to calculate the scattering coefficient u al 
and the transport scattering coefficient u s1 ' of the scattering medium. This sort of the arithmetic processing 
is operated at high speed by a microcomputer installed in the arithmetic processing unit 16 

In the above case, a light source, a light incident means, a photodetecting means, an amplifying means 
for a detected optical signal which are described in the first embodiment can be utilized. These apparatus are 
to measure the optical constants of an unknown scattering medium and for general purpose, which can be util- 
ized in various fields such as science and engineering, biology, spectral analysis, and others. 

(2.3) Third Embodiment 

The third embodiment of the apparatus carrying out the method of measuring internal information in a scat- 
tering medium according to the present invention is shown in Fig. 21. This apparatus is for the measurement 
and/or the mon.tonng of the concentration of oxyhemoglobin or the degree of oxygen saturation of hemoglobin 
in a human brain. A band-type container 26 is put on a head 30 like a headband to measure the degree of 
saturation of hemoglobin in a brain. The apparatus of this embodiment uses light having two predetermined 
wavelengths X, and X 2 , and the operation is substantially the same as the first embodiment In Fig 21 the 
parts each which has the same function as described in Fig. 7 are represented by the same reference numerals 

A wavelength of Lights having two different wavelengths X, and X 2 are selected by a wavelength selector 
4. and light is incident on the head 30 which is an object to be measured through a light guide 6 The two wa- 
velengths should properly be selected. Light diffused during propagation in the head which is an object to be 
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measured is received by light guides 8, 9 and 10 arranged at positions spaced apart from the light incidence 
position by distances p 1f p 2 , and p 3 , and converted and amplified into electric signals by the first photodetector 
12, the second photodetector 13 and the third photodetector 14. In this case, a power supply, or a signal is 
coupled to external equipments (not shown) through a connector 28 installed in a container 26. The arithmetic 

5 processing (not shown) to signals is the same as the first embodiment, and in this embodiment, the degree of 
oxygen saturation of hemoglobin is displayed. In addition, the signals may be converted into electric wave or 
optical signals and then transmitted. 

A light source, a light incidence means, a photodetecting means, an amplifying means for the detected 
optical signals or others which are described in the first embodiment can be utilized. Surface reflection or a 

10 space between the light guide and the head may be a problem for the human head. In this case, the interface 
material may be utilized as described above. Fig. 22 shows a modified example of a light incident position and 
a detector shown in Fig. 21 in detail. In this example, light guides shown in Fig. 21 are omitted, and the interface 
material 11 having substantially the same scattering coefficient and the absorption coefficient as an object to 
be measured is used in a space between a head and a wavelength selector 4, and a space between a head 

15 and a photodetector 12. In this way, a difference of the states of light coupling between photodetection at p 1( 
p 2 and p 3 and the head can be reduced, and efficiency of light collection is improved. 

These apparatus can be used, e.g., for the measurement or the monitoring of the concentration of oxyhe- 
moglobin in a muscle of a human who runs a marathon. 

As described above, according to the method of measuring internal information in the scattering medium 

20 and the apparatus for the same, the measurement of the internal information in the scattering medium at high 
precision, that is, the absolute values of the transport scattering coefficient u s * and the absorption coefficient 
u, a can be measured. If the principle of the dual-wavelength spectroscopy is further utilized, the internal infor- 
mation such as the concentration of the specified constituent can accurately be measured. The spacial dis- 
tribution and time-change in the internal information can be measured. Moreover, since the measurement ap- 

25 paratus which utilizes the present invention uses the time-integrated value of the optical signal, light utilization 
factor is high, and as the signal-to-noise ration is large, the measurement precision becomes high. Therefore, 
the non-invasive measurement and imaging for the internal information of a human head, a body, a plant such 
as a standing tree, or others can be operated by a simple apparatus. 

30 

Claims 

1. A method of measuring a scattering medium comprising: 

a step of causing light having a predetermined wavelength to be incident on the scattering medium; 
35 a step of detecting light having said predetermined wavelength diffused during propagation in said 

scattering medium at detection points corresponding to three or more kinds of points where the distance 
between a light incidence position and a photodetection point is different to obtain the three or more kinds 
of the detected signals; and 

a step of extracting internal information in said scattering medium by performing arithmetic proc- 
40 essing to said detected signals based on three or more kinds of simultaneous relations between a scat- 

tering characteristic and an absorption characteristic on a diffusion-propagation path in said scattering 
medium when the light having said predetermined wavelength is diffused during propagation and reaches 
at said detection points, and said detected signals. 

45 2. A method of measuring a scattering medium according to Claim 1, wherein light to be incident on said 
scattering medium is pulsed light; and said simultaneous relation is a following equation, 

lnQ[pJ = ln[TiZo/(2tt)] - (3/2)^* + pfl + In^Zo 2 + Pi 2 ) 1 * + 1] - Hefffo 2 + Pi 2 ) 1/2 
for 



Zo = a% a + us') 

so where 

Q: the measured light intensity 

p,: the distance between the light incidence position and the photodetection point 

i: the integer from one to an integer above three 

q: the light attenuation factor when light is incident on the scattering medium and when light is 

55 emerged 

u^ffC the effective attenuation coefficient 

(m*f = [3^ (Pa + us')l 1/2 
a: the constant determined in accordance with a kind of the scattering medium 
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P*: the absorption coefficient 

Ps*: the transport scattering coefficient 

3. A method of measuring a scattering medium according to Claim 1, wherein said light having said prede- 
termined wavelength is light having two or more kinds of wavelengths, each having a different absorption 
coefficient to a specified absorptive constituent of said scattering medium; and 

said arithmetic processing comprises a process of extracting primary internal information in said 
scattering medium obtained by the measurement with one kind of said predetermined wavelength, and a 
process of extracting secondary internal information in said specified constituent of said scattering me- 
dium, utilizing said primary internal information corresponding to light having said two or more kinds of 
wavelengths. 

4. A method of measuring a scattering medium according to Claim 1 , wherein said internal information is an 
absorption coefficient and a transport scattering coefficient. 

5. An apparatus for measuring a scattering medium comprising: 

light-emitting means for emitting light having a predetermined wavelength; 
light-incident means for causing the light having said predetermined wavelength to be incident on 
the scattering medium; 

photodetecting means for detecting light having said predetermined wavelength, diffused during 
propagation in said scattering medium at photodetection points corresponding to three or more kinds of 
points where the distance between a light incidence position and a photodetection point is different, and 
obtaining three or more kinds of the detected signals; and 

arithmetic processing means for extracting internal information in said scattering medium by per- 
forming arithmetic processing to said detected signals, based on three or more kinds of simultaneous re- 
lations among a scattering characteristic and an absorption characteristic on a diffusion-propagation path 
in said scattering medium when the light having said predetermined wavelength is diffused during prop- 
agation and reaches at said detection points, and said detected signals. 

6. An apparatus for measuring a scattering medium according to Claim 5, wherein said light-emitting means 
emits pulsed light; and said simultaneous relation is a following equation, 

lnQ[pJ = ln[TiZo/(2*)] - (3/2)ln(zo 2 + Pi 2 ) + \nM*> 2 + Pi 2 ) 172 * 1] - p^Zo 2 + Pi 2 ) 172 
for 

zo = a/(u. a + n s ') 

where 

Q: the measured light intensity 

P S : the distance between the light incidence position and the photodetection point 

i: the integer from one to an integer above three 

n: the light attenuation factor when light is incident on the scattering medium and when light is 

emerged 

Meff: the effective attenuation coefficient 

<Peff = [3p a (u a + us')]"* 
the constant determined in accordance with a kind of the scattering medium 
m*: the absorption coefficient 

u. s *- the transport scattering coefficient. 

An apparatus for measuring a scattering medium according to Claim 5, wherein said light having said pre- 
determined wavelength is light having two or more kinds of wavelengths, each having a different absorp- 
tion coefficient to a specified absorptive constituent of said scattering medium; and 

said arithmetic processing means executes a process of extracting primary internal information in 
said scattering medium obtained by the measurement with one kind of said predetermined wavelength 
and a process of extracting secondary internal information in said specified constituent of said scattering 
medium, utilizing said primary internal information corresponding to light having said two or more kinds 
of wavelengths. 

A method of measuring a scattering medium according to Claim 5, said internal information is an absorp- 
tion coefficient and a transport scattering coefficient. 

A method of measuring a characteristic of a sample comprising irradiating the sample at a position of in- 
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cidence and detecting the radiation after passage through the sample at three detection positions, each 
detection position being at a different respective distance from the position of incidence. 

10. A method as claimed in claim 9, wherein the characteristic comprises a scattering characteristic and/or 
an absorption characteristic, the method further comprising generating signals in response to the radiation 
detected at the three detection positions and processing said signals thereby to obtain a measure of said 
characteristic. _ 
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